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n INTRODUCTION 


The purpose of this appendix is to derive expressions for computing the Aerody- 
namic influence coefficient matrix (or [A] matrix) for nonplanar wing-body-tail configura- 
tions. An aerodynamic influence coefficient ajj (i.e. an element of the A matrix) is 
defined as the load in lbs. (Newtons) induced on panel i as a result of a unit angle of 
attack (rad) on panel j. Simulation of wing-body interference is not attempted . Fuselage, 
wing and tail thickness are assumed to be small with the result that the thickness effect on 
the flow-field is negligible. No dihedral effects are considered, even though these can 
be easily Included by computing both downwash and sidewash. Camber effects may also 
be included. Symbols used in this appendix are defined in Section 2. 

The method for determining the aerodynamic influence coefficient matrix is 
based on the "lifting" solution to the small perturbation, steady potential flow equation. 
Expressions for the various velocity functions, needed for computing the downwash are 
presented in Section 3 of this appendix. In section 4, a comparison of pressure distribu- 
tions and stability derivatives derived from this investigation, with already existing results 
is made . 



2. SYMBOLS 


The units used for the physic.al quantities defined in this paper are given both In 
the International System of Units (SI) and the U.S. Customary Units. 


Symbol 

Description 

Dimension 

a 

Downwash influence coefficient 

Nondime nsional 

[A] 

Aerodynamic influence coefficient 

ft^ rad ^ (m^ rad 

AR 

Aspect ratio 

Nondimensional 

bjk 

L j k //3 

Nondimensional 

Cp 

Pressure coefficient 

Nondimensional 

Fa 

' AerodynamTc Pressure force on a panel 

Lbs. (Newton) 

F 

Velocity function 

Nondimensional 

K 

Constant 


L 

Slope of leading or trailing edges of a 
panel 

Nondimensional 

Moo 

Free stream Mach number 

Nondimensional 

q 

Dynamic pressure 

Lbs. ft ^ (Newton m ^) 

s 

Area of a panel 

ft^ (m^) 

U 

Perturbation velocity in x-direction 

Nondimensional 

V 

Perturbation velocity in y-direction 

Nondimensional 

w 

Perturbation velocity in z-direction, 
downwash 

Nondimensional 

(X/ // z) 

Rectangular Cartesian coordinate system 

ft. (m) 

z 

c 

Ordinates of Camber line 

ft. (m) 


2 





rad 

rad 



Nondimensional 


( v) 


Velocity potential Nondimensional 

Integration variables in Cartesian system ft, (m) 



Subscripts 

i 

j 

k 

k 

( )■ 

Matrices 

[ ] 

1 ! 

1 f 

N 


Spanwise station in fraction of the semi-span 

Aerodynamic panel number 

Aerodynamic panel number 

Panel corner point 

Aerodynamic panel number 

Referred to primed system of coordinates 

Square matrix (n x n) 

Column matrix (n x 1) 

Row matrix (1 x n) 

Diagonal matrix (n x n) 
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3. METHOD OF ANALYSIS 


The aerodynamic influence coefficient method is based on the "lifting" solution 
to the small perturbation, steady potential flow equation; 


(1 - + <l> = 0 

V yy zz ( 1 ) 

where 0 is the perturbation velocity potential and the freestream Mach number. 

Note that the same equation is also satisfied by the perturbed velocity components. It is 
well known that the lifting solution of the above equation can be obtained with a vortex 
distribution with strength to be determined by satisfying a boundary condition of flow 
tangency on the wing surface. In this appendix the Woodward's method of solution 
(Reference 1), has been used. 

In Woodward's method, the entire wing or airplane planform is divided into a num- 
ber of quadrilateral panels os shown in Figure (1). On each panel, the Incidence and 
camber effects are represented by a constant planar pressure vortex distribution with an 
unknown strength A u. The velocity potential at any point in the flow field due to this 
vortex distribution on one panel is given by: 


</» (x', y'. 



Panel 


z'(x' -^) d^dT, 

[( y ' - [(y. . ,)2 ^,2 , 


^ (1 - K)An 

2 7T 


f r z'd^dii 

jj (y* - 7,)^ + z'2 


( 2 ) 


where K ~ 0.5 for M^< 1 and K = 1.0 for Mj^ > 1. The primed coordinate system is 
a local coordinate system as shown in Figure (2), There are many ways of finding the in- 
tegrals in Equation (2). A simple way is to regard each panel as an algebraic sum of four 
semi-infinite triangular regions as shown in Figure (3a). 

Thus, the velocity potential at any point P Is given by: 

0 (P) = 0-| (semi-infinite triangular region from corner 1) 

- <^2 (semi-infinite triangular region from corner 2) 

/ 

~ 4> 2 (semi-infinite triangular region from corner 3) 

(semi-infinite triangular region from corner 4) 
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W' 



Figure 1 Example of Quadrilateral Panel Distribution 
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A similar approach can also be used when the leading and trailing edges of the 
panel have negative slopes and when points 2 and 4 are coincident, i.e., the panel is 
triangular. In the latter case, the same scheme as for the positive leading edge shape can 
be used. However, when the leading edge slope is negative, the scheme shown in Figure 
(3b) is convenient. 

It is clear then that the fundamental region of integration for Equation (2) is the 
semi-infinite triangular region with a positive leading edge slope. If the equation of the 
leading edge is given by ^ = Lrj , the velocity potential due to a constant pressure vor- 
tex distribution in the semi-infinite triangular region from the corner k is as follows: 


A 


KA 


^I 

f 


z:' (x* - d^ dr? 


0 - K)Au 

21r 


[ (y- - (x- - ^)2 + y. > ^)2 + ^.2^ 




z'd^dT? 

Lt 7 (y' -r}) + z'^ 


(3) 


where x* = - X|^, y' ^ Xp ■" /[<. ~ In subsonic flow, = oo and i?] 

In supersonic flow, the integration limits are as shown in Figure (4) and are given below: 

<2 , \ r , . .2 


00 


“ {(^oo".^) i(y' - vr + z'^] } 


(3a) 




Lx' - (M^ - 1) ^ y' - { (M^- 1) [(x' - Ly')^ + (L^ + 1 - z'^ ] 

+ 1 - 


(3b) 
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Equations (3) have been integrated In closed form by Woodward (Reference 1) 
The result Is: 


^ _ KAu 

*k ■ "27* 


{ (x' - Ly') (Fj + F^) + z' 


F 

[r+l-M^ -|-L(F, -Fj) ij 


(4) 


The corresponding velocity components are: 


/ X _ _ KAu (r: \ 

“k <P) - a; 57 <^3 + ^ 4 > 

P 

/ X_ ^^k _ KAu r i /c 4 -c \ ,oc 1 

'^k " IT " ~ ^ ^ 

'P 


(5) 


/x_^'^k_KAu 

p 


r -,2 


+ 1 - M: 


00 


F,-L(F, - FJ-y'SF, 


/3 z ■ 1 '5' 

The velocity functions F's are defined as follows: 


=Re { log 




x' + d' 


1 - r' 


F, = /3 Re I log v - } 




2 2 
r + 1 - ) 


Fo = Re I tan ^ 

. Ir'^^-xy 


Fz = tan’"^ ^ 
4 z' 


for < 1 


F4=G 


for Mqq > 1 


(6a) 


(6b) 


(6c) 


(6d) 
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o> 

11 

LL. 

forM <1 

CO 

r 

o 

11 

LL. 

for >1 

^ x'+d' 

for 

F - 

^6 y3r'2 

for Mqq>1 


2 2 '/2 

where: r* = (z' + y‘ ) 

2 2 2 1 
d' = tx‘ + (1 - M^) r'^ 1 

T = Lx' + (1 - Mi)r 

'y~ — x‘ - Ly‘ 

_ 2 2 ' 2 2 
r = I(x' - Ly')^ + (L"^ + 1 - M^)z'^ ] 

d =l3r + 0 - M^r ] 


(6e) 


(6f) 

(7a) 

(7b) 

(7c) 

(7d) 

(7e) 

(7f) 

(7g) 


Let p be the downwash control point of a panel . Then the downwash at p of panel 
i due to a unit vortex distribution on the semi-infinite triangular region from the corner k 
of panel j is given by: 

"ijk W ' ^ i i - '"ik ‘f 1 ■ >'ijk 


where the plus sign is for the subsonic case and the minus sign for the supersonic case. The 

quantities b* and y' in Equation (8) are defined as follows: 
jk ' ijk 
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for k = 1 / 2 


b- = ilL = i. 

jk p /3 yj2 -yji 

= 1 

^ yj4 - y j3 


for k = 3, 4 I 


% = yryjk 


( 9 ) 


00 ) 


x'..|, and z'jjb defined in a similar way as; 

11 

_x 

1 


z'..i = Z. - Z.| 
ijk 1 jk 

( 12 ) 

and also 

11 

(13) 


Equafion 8 Indicates that only , ? 2 t ^^5 velocity functions are needed 

to evaluate the downwash W|^ . Simplified expressions of these velocity functions can be 
derived by using equations ^7 and 9 through 13. These expressions for both subsonic and 
supersonic cases are presented in Table 1 . 

Referring back to equation 8 , it follows that the downwash at point p of panel i 
due to a unit vortex distribution on the jth panel, represented by a .., , can be written 
as: 


a 


ij1 






( 18 ) 


In this report oni/ symmetric flight configurations are considered. Let the mirror image 
of the jth panel with respect to the line of symmetry be denoted by j'^'^ panel . The down- 
wash at point p due to the panel j' is aj* 2 . The total downwash at p due to a unit vortex 
distribution on both j<^h and panels ’Is then: 


a.. 


U 



(19) 
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Table 1 Velocity functions for computing downwash 
according to equation 8 


Subsonic Case; 

FI = Reln 






|8 




Supersonic Case: 



(14d) 


13 





F5 = 0 



F2 = 0 



ify' < b'^» 
■fy' = b'g' 


(15e) 

(15f) 


( 16a) 
( 16b) 
( 16c) 



F5 = 0 


Ify' > b'^‘ 


F6 = 0 


Nofe: If = yy^ + z'^ and b' Vl, all functions are zero. 


(16d) 

(I6e) 

(16f) 



Case 3. 




and b'< 1 . 


FI = 0. 


F2 = 0. 


F2 = 0 


F2 = 


7T 

2\l\ - b‘^ 


F2 = 





(17a) 

ify' 0 


(17b) 

if y' )>• 0 and < (by 

1 - b'^lz'l ) 

(17c) 

if y' y 0 and = (b‘y' +\/ 

> 2 

1 -b'^(z'l) 

(17d) 



(17e) 


if y' > 0, ) (by +p - b'^ |z‘ I ) and y‘> b'^' 


F2 = 


ij] - b '^ 


F2 = 0 


07f) 

if y' > 0, ^ ^ (b‘y‘ "V ^ “ b'^ j z'j ) and y' = b’ 

( 17g) 

if y' > 0, r' > (by +\A - b'^[ z' 1 ) and y'< b' | ' 


F5 = 0. 


(17b) 


Note: 


F6 = 0. 

If ?'< \/y'^ + z'^ and b' ^ 1 , all functions are zerOo 


( 17i) 
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The boundary conditions require that the flow be tangential to the camber line. 
Since the downwash at p of panel i due to the vortex distribution on all panels is: 


N 


w 


E a.. 

ij 

j=1 


(Au), 


where Ny^ is the number of panels on one^half of the wing, fuselage and tail , the tan- 
gency condition can be written as: 


N 


dz 

( — 
dx 


w 


- a )| - °ij 


(20) 


or in matrix form: 


dz 

1 " 


ax 


a I = [a.. ){Au| 


( 21 ) 


Solving this equation for |Au| , the unknown vortex strengths Au are obtained as: 

( 22 ) 


{Au| = [a., j - oc\ 

dx 


Note that Au = u"^-u . Since C = -2u, it follows that AC - G -C - 2Au . 

P P P P 

Hence: 

|ACp! = 2 la,.]-' 1-^ - (23) 

The total force on any panel i is given by; 



where S; is the area of panel i andq is the' dynamic pressure. In matrix form, the panel 
forces can be written as: 
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0 


{FaH 




0 • s 


N 




w 


= 2q 


N1 

^ [«w} 


.,r' 


( 25 ) 


Where [A] 


is the aerodynamic influence coefficient matrix and is given by: 


J 



( 26 ) 


and 


a 


w 


dz 

c 

dx 


( 27 ) 


The computer program for generating this [A ] matrix is included in Reference 2 
(Appendix A of the Summary Report). 
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4. VALIDATION OF THE [A] - MATRIX ROUTINE 


To check out the computer program for the aerodynamic matrix, equation (23) was 
used to compute the pressure distribution for a highwing-midtail configuration, shown in 
Figure 5. R^ults were compared with those derived from a NASA Langley Vortex Lattice 
program. Reference 3. 

Both programs employed a total of 80 panels with 50 panels on the halfwing, 10 
on the half fuselage and 20 on one side of the tail . The pressure distributions are com- 
pared in Figure 6. Note that the computed AC„ is plotted at the panel centroid for the 
K.U. program, but at the panel quarter chord for^ the Vortex Lattice program. It is seen 
that the results from both programs agree reasonably well , The computation of stability 
derivatives is compared in Table 2. The experimental data are obtained from Reference 4. 
To determine the effect of the number of panels on the numerical outcome, 100 and 120 
panel solutions were computed. Table 2. shows that there was little change by using a 
panel scheme of more than 100. 


Table 2. Comparison of Computed 
Stability Derivatives with Experiment for 


the Configuration of Figure 5 at M - 0.25 



' 

K.U. Program 

Vortex Lattice Program 

Experiment 
• Reference 4 

Derivatives 

C. rad ^ 
•-« 

No. of panels 

No. of panels 

80 

100 

120 

80 

3.28 

3.1930 

3.105 

3.0780 

3.1208 


-0.0395 

-0.0541 

-0.0429 

-0.0438 

-0.065 


In comparison with experimental results, it is assumed that for a high wing con- 
figuration, the wing may be placed directly above the fuselage at a distance of one half 
fuselage diameter, so that two lifting surfaces (wing and fuselage) will overlap in some 
region. A similar aerodynamic representation has been assumed for a low wing configur- 
ation or for a high or low tail . Several other configurations presented in References 5 
through 8 were investigated with the K.U. program. Their geometries are shown in Figures 
7 through 10 respectively. The total number of panels used on a half configuration was 
100. Table 3, shows the comparison with experimental results. It is seen that the com- 
puted results show reasonable agreement with experiments. 
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WING ROOT CHORD 
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Figure 5 Geometry of a. wing-body-tail configuration used for program checkout 




















ALL BRACKETED DIMENSIONS IN METERS 

NOTE; CONFIGURATION WBH 2 OF 
REFERENCE 5 


Figure 7 Aerodynamic RepresentaHon for Geometry of WBH2 
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Figure 8 Aerodynamic Representation for Geometry of WBH 3 




c/4 LINE 


ALL BRACKETED DIMENSIONS IN METERS 
NOTE; CONFIGURATION WBH 4 OF REFERENCE 7 


Figure 9 Aerodynamic Representation for Geometry of WBH 4 
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NOTE' CONFIGURATION WBH 5 OF REFERENCE 8 ‘ (0.0443) 



Figure 10 AerodynamiG Representation for Geometry of WBH 5 











Table 3 Comparison of Computed Aerodynamic Characteristics with Experimeni'al Results 
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